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ABSTRACT 



We introduce a self-similar mass model for early-type galaxies, and constrain it using 
the aperture mass-radius relations determined from the geometries of 22 gravitational 
lenses. The model consists of two components: a concentrated component which traces 
the light distribution, and a more extended power-law component (p oc r~ n ) which 
represents the dark matter. We find that lens galaxies have total mass profiles which 
are nearly isothermal, or slightly steeper, on the several-kiloparsec radial scale spanned 
by the lensed images. In the limit of a single-component, power-law radial profile, 
the model implies n = 2.07 ± 0.13, consistent with isothermal (n = 2). Models in 
which mass traces light are excluded at > 99% confidence. An n = 1 cusp (such as 
the Navarro- Frenk- White profile) requires a projected dark matter mass fraction of 
fcdm = 0.22 ± 0.10 inside 2 effective radii. These are the best statistical constraints yet 
obtained on the mass profiles of lenses, and provide clear evidence for a small but non- 
zero dark matter mass fraction in the inner regions of early-type galaxies. In addition, 
we derive the first strong lensing constraint on the relation between stellar mass-to-light 
ratio T and galaxy luminosity L: T oc L ' -°- 12 , which is consistent with the relation 
suggested by the fundamental plane. Finally, we apply our self-similar mass models to 
current problems regarding the interpretation of time delays and flux ratio anomalies 
in gravitational lens systems. 

Subject headings: galaxies: elliptical and lenticular - galaxies: structure - gravitational 
lensing 



1. Introduction 



The structure of galaxies is closely linked to their formation and evolution, and therefore 
provides a vital testing ground for the cold dark matter (CDM) paradigm which has proven so 
successful on large scales (e.g., Spergel et al. 2003). The primary observable is the shape of the 
radial mass profile, from which the relative distributions of luminous and dark matter can be 
determined. The outermost regions of galaxies have now been extensively probed using weak lensing 
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(e.g., Fischer et al. 2000; McKay et al. 2001; Kleinheinrich et al. 2003) and satellite dynamics (e.g., 
Zaritsky et al. 1997; Romanowsky & Kochanek 2001; McKay et al. 2002; Prada et al. 2003), each 
of which reveals an extended dark matter halo out to 100 — 300 kpc. Methods for tracing the inner 
regions of galaxies largely depend on the galaxy morphology. Spiral galaxies are the most easily 
studied, and are consequently the best understood, as rotation curves can be mapped to ~ 30 kpc 
using kinematic tracers such as HI. Flat rotation curves at > 5 kpc clearly indicate the presence of 
dark matter (e.g., Rubin, Thonnard & Ford 1980), but dynamical observations also suggest that 
the inner halo may be significantly less concentrated (e.g., McGaugh & de Blok 1998; de Blok et al. 
2001; Salucci 2001; de Blok & Bosma 2002) than the cuspy halos predicted by CDM simulations 
(e.g., Moore et al. 1999b; Bullock et al. 2001). 

Compared to their late-type counterparts, early-type (elliptical and lenticular; E/S0) galaxies 
are not as well understood. Some of the best evidence for dark matter outside of a few optical radii 
has been obtained from investigations of the extended, X-ray emitting gas (e.g., Fabbiano 1989; 
Matsushita et al. 1998; Loewenstein & White 1999). The inner regions of early-type galaxies are 
typically probed using stellar dynamics (e.g., Rix et al. 1997; Gerhard et al. 2001). However, the 
velocity dispersion is not an unambiguous tracer of the mass profile, as it requires an understanding 
of the orbital anisotropy. In addition, the existence of a tight fundamental plane (FP; Djorgovski 
k, Davis 1987; Dressier et al. 1987) presumably provides clues about the inner structure of early- 
type galaxies, but little consensus has been achieved regarding its meaning (e.g., Faber et al. 1987; 
Renzini & Ciotti 1993; van Albada, Bertin & Stiavelli 1995; Ciotti, Lanzoni & Renzini 1996; Pahre, 
de Carvalho, & Djorgovski 1998; Bertin, Ciotti, & Del Principe 2002). Many outstanding questions 
therefore remain regarding the structure of E/S0 galaxies, the universality of this structure, and 
its relation to the predictions of CDM models. 

Strong gravitational lensing is a powerful and unique tool for investigating the mass distri- 
butions in early-type galaxies at intermediate redshift, as it probes mass directly. First, model- 
independent projected masses may be inferred from the geometry of the lensed images (e.g., Schnei- 
der, Ehlers & Falco 1992). When mass models are normalized by this constraint, the stellar velocity 
dispersion becomes a far more sensitive probe of the radial mass profile. Ongoing efforts by the 
Lens Structure and Dynamics (LSD) survey (Koopmans & Treu 2002, 2003; Treu & Koopmans 
2002a; collectively "KT") have made significant progress in obtaining and interpreting dynami- 
cal measurements in gravitational lens systems. Second, models of the galaxy mass distribution 
can be constructed based on the positions, morphologies and flux ratios of the lensed images. Of 
particular importance are the additional astrometric constraints obtained from complex structure 
in the lensed source, which can be used to break degeneracies and distinguish among radial mass 
profiles. These constraints include large-scale extended emission from radio lobes (e.g., Kochanek 
1995) or quasar host galaxies (see, e.g., Kochanek, Keeton, &: McLeod 2001 and Saha & Williams 
2001 for general discussions), lensed images of multi-component sources (e.g., Cohn et al. 2001; 
Muhoz, Kochanek & Keeton 2001), and the relative orientations of milliarcsecond-scale radio jets 
(e.g., Rusin et al. 2002; Winn, Rusin & Kochanek 2003). 
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The mass distributions in several lens galaxies have now been investigated using either direct 
modeling or stellar dynamics. In all but one case, the profile is consistent with an isothermal 
model (p oc r~ 2 ). The lone exception is PG1115+080, for which the high stellar velocity dispersion 
(Tonry 1998) suggests a mass profile steeper than isothermal (p oc r~ 2 - 35 ; Treu & Koopmans 
2002b), although this also implies that the galaxy lies well off the FP. The lensing results are 
therefore generally consistent with the nearly isothermal profiles favored by X-ray and dynamical 
studies of local elliptical galaxies. 

Understanding the mass profiles of lens galaxies is required for the determination of the Hubble 
constant (Hq = 100/ikms" 1 Mpc -1 ) from measured time delays. Kochanek (2002) demonstrates 
that the predicted delay between two images, and hence the derived value of Hq, is primarily 
governed by the average surface mass density in the annulus defined by their radial distances from 
the lens center, with a small correction for the slope of the profile in this annulus. Each of these 
quantities is determined by the shape of the radial mass profile. Based on a study of several lens 
systems which are dominated by a single galaxy and have well-measured time delays, Kochanek 
(2002) finds h = 0.51 ± 0.05 if lenses are isothermal. This value is lower than those obtained 
from the Hubble Space Telescope (HST) Key Project (h = 0.72 ± 0.08; Freedman et al. 2001) and 
WMAP data (h = 0.72 ±0.05, although this requires the assumption of a flat cosmology and a dark 
energy equation of state of w = —1; Spergel et al. 2003). In fact, Kochanek (2002, 2003) shows 
that early-type galaxies must be described by models in which mass traces light, and hence have 
no extended dark matter halos, to reconcile current time delay measurements with h ~ 0.7 (see 
also Williams & Saha 2000). This apparent discrepancy with other observational constraints and 
the CDM paradigm demands continued examination of the lens galaxy population. 

Decomposing the luminous and dark matter components of early-type galaxies is also poten- 
tially important for investigating the nature of flux ratio discrepancies in gravitational lens systems. 
It has now been extensively demonstrated that simple, smooth mass models have great difficulty 
reproducing the observed flux ratios in four-image lenses (e.g., Dalai & Kochanek 2002; Metcalf 
& Zhao 2002; Chiba 2002; Keeton, Gaudi & Petters 2002), both in the optical and radio. The 
effect has been traced to small-scale structure in the gravitational potential (Kochanek & Dalai 
2003), and can be due to either stars (microlensing; e.g., Witt, Mao, & Schechter 1995; Schechter 
Sz Wambsganss 2002) or CDM satellite halos (millilensing; e.g., Mao & Schneider 1998; Metcalf 
& Madau 2001; Bradac et al. 2002; Chiba 2002; Dalai k Kochanek 2002). Both phenomena arc 
surely present, as radio sources are generally too large to be affected by stars, while uncorrelated 
time variability is an unmistakable signature of stellar microlensing in some optical lenses (e.g., 
Wozniak et al. 2000). Microlensing is most efficient at accounting for observed anomalies, in partic- 
ular the preferred demagnification of the saddle point image, when stars comprise a small fraction 
(15-30%) of the surface mass density at the image radius (Schechter &: Wambsganss 2002). If we 
can constrain a two-component model of the galaxy mass distribution, then we can predict the 
stellar surface mass fraction at the radii typically spanned by lensed images. Thus, the predictions 
and requirements of the microlensing hypothesis may be tested. 
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Because mass distributions have been directly constrained in only a small number of gravita- 
tional lenses, it is worthwhile to consider statistical constraints on the lens galaxy population. Such 
attempts have thus far been very limited. Several analyses (e.g., Rusin & Ma 2001; Evans & Hunter 
2002) have focused on the absence of central images in deep radio maps of gravitational lens sys- 
tems, which has been invoked to rule out profiles that are much shallower than isothermal or have 
significant cores. Keeton (2003) suggests that such constraints have been over-interpreted, since 
the magnification of the missing image is strongly dependent on the properties of the mass profile 
very close to the galaxy center. In addition, because mass distributions steeper than isothermal do 
not produce central images, any constraints based on the absence of these images are one-sided. 
The number and sizes of lens systems have also been used to investigate galaxy mass profiles in the 
context of CDM models (Keeton 2001). It is shown that adiabatic compression (e.g., Blumenthal 
et al. 1986) yields nearly flat rotation curves for a range of initial halo concentrations and cooled 
baryon fractions (see also Kochanek & White 2001 and Kochanek 2003), but typical CDM halos 
(e.g., Bullock et al. 2001) appear to be too centrally concentrated to account for the data. The 
analysis, however, offers few quantitative results describing the mean mass profile and scatter, and 
therefore demands further study of two-component mass models. 

In this paper we introduce a new technique to statistically constrain the mass profiles of early- 
type lens galaxies via aperture mass measurements. Each lens must obey a strict relationship 
between the image radii and the projected masses they enclose. The set of mass-radius relations 
can be used to investigate the typical mass profile in a sample of lenses if we assume that early-type 
galaxies have a self-similar (homologous) structure. This means that the functional form (shape) of 
the mass distribution is the same from galaxy to galaxy, and is scaled by the properties of the light 
distribution. An underlying homology between mass and light is suggested by the existence of a 
tight fundamental plane relating optical and dynamical observables. In one popular interpretation 
of the FP (e.g., Faber et al. 1987; van Albada et al. 1995), early- type galaxies are structurally 
homologous, but the mass-to-light ratio has a luminosity dependence. Gerhard et al. (2001) use 
dynamical data to bolster this hypothesis, demonstrating that more luminous early-type galaxies 
have higher stellar mass-to-light ratios in their central regions. Other analyses have cast doubt on 
the description of the FP in terms of strong homology and varying mass-to-light ratios (e.g., Caon, 
Capaccioli, Sz D'Onofrio 1993; Pahre et al. 1998; Bertin et al. 1994, 2002), and instead suggest 
a weak homology in which structural parameters vary systematically with luminosity. However, 
because a strong homology is the simplest reasonable assumption, it is a good starting point for 
the statistical investigation of mass profiles in early-type galaxies. 

Section 2 details the aperture mass-radius constraints offered by lensing. Section 3 describes 
the lens sample and outlines our various assumptions. Section 4 introduces our self-similar models 
and the associated methodology. Section 5 presents our statistical constraints on the radial mass 
distribution in lenses. Section 6 demonstrates how these constraints can be used for the interpreta- 
tion of measured time delays and flux ratio anomalies among lensed images. Section 7 discusses our 
findings in the context of recent astrophysical and cosmological results. Unless otherwise stated, 
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we assume a flat VLm = 0.3 cosmology with h = 0.65 for all calculations. 



2. Aperture Masses from Lens Geometries 



The geometries of gravitational lens systems yield model-independent constraints on projected 
masses. Consider, for example, a deflector with a circularly symmetric surface density distribution, 
and a source which sits directly behind it. The image will form on a ring of (physical) Einstein 
radius REin, which is related to the aperture mass [M-Ein) it encloses by 



where S cr is the critical surface density. Angular diameter distances to the lens, to the source, 
and from the lens to the source are Dd, D s and Dd s , respectively. For this simple case, the mass 
enclosed by the ring is determined by its radius alone, and is therefore independent of the shape of 
the mass profile. 

Real lenses are more complicated. First, sources (and their images) are, in general, not sym- 
metrically placed with respect to the mass distribution. This means that the geometry may be 
described by more than one characteristic radius. Second, gravitational potentials are, in general, 
not circularly symmetric. The presence of a quadrupole, due to either external shear or ellipticity, 
effectively stretches the image plane along one axis and compresses it along a perpendicular axis, 
thereby confusing the interpretation of geometry in terms of mass. These issues are surmountable, 
however, allowing any lens geometry to yield model-independent relations between masses and 
radii. 

Four-image lenses (quads) offer the best mass constraints, as the above complications are 
minimized. The images typically reside at similar radii, and their good angular coverage allows 
for the robust determination and removal of the quadrupole. Just about any mass profile can 
reproduce the image positions and fluxes in a quad, but modeling demonstrates that there is a 
radius at which all models agree on the enclosed projected mass (e.g., Cohn et al. 2001; Muhoz 
et al. 2001). We have tested the generality of this claim using Monte Carlo simulations. First, 
we create fake four-image lenses from an elliptical mass distribution with some fixed radial profile 
residing in an external shear field of random orientation. We then model the image positions and 
fluxes using a broad range of mass profiles, tabulating the mass enclosed by a circular aperture as 
a function of radius for each best-fit model. We find that the suite of aperture mass curves M(R) 
cross at nearly the same radius (REin), the value of which can be accurately estimated (to ~ 1%) 
from the critical radius of a model-fit singular isothermal sphere (SIS) in an external shear field. 
This is true regardless of whether the mass distribution is spherical or elliptical. The enclosed, 
projected mass inside REin is then given by eq. (1), to an accuracy of ~ 2%. 

Two-image lenses (doubles) yield weaker mass constraints. If the potential is circularly sym- 
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metric, then the lens equation dictates a model-independent mass-radius relation: 



1 [ M(fli) M(R 2 ) 
tt [ Ri + R 2 
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where i?i and i?2 are the radii of the images with respect to the galaxy center, and M(R\) and 
M(i?2) are the projected masses they enclose. Hence, doubles constrain a combination of two radii 
and two masses. Note that in the limit of a ring {R\ = R2), we recover eq. (1). The presence of 
a quadrupole smears the mass-radius relation for doubles, and the limited number of constraints 
means that its effects cannot be removed in a model-independent manner. If the quadrupole is 
represented as an external shear, then the maximum fractional deviation from eq. (2) is equal to 
the shear magnitude 7, and this occurs when the shear axis is parallel or perpendicular to the axis 
defined by the lensed images. Averaging over orientations, the rms deviation is ~ 7/V2. We have 
confirmed this assertion using Monte Carlo simulations, in which fake doubles are produced from 
one mass profile and modeled using a range of mass profiles. While four-image lenses typically 
have quadrupoles which can be described by shear fields of amplitude 0.1 < 7 < 0.2 (e.g., Keeton, 
Kochanek & Seljak 1997; Holder & Schechter 2003), we might expect doubles to have somewhat 
lower quadrupoles, on average, as they represent a greater fraction of the lensing cross section 
when the shear is small (e.g., Keeton et al. 1997; Rusin & Tegmark 2001; Finch et al. 2002). 
Consequently, deviations from eq. (2) should generally be ~ 10%. We therefore employ eq. (2) 
as a model-independent mass-radius relation for doubles, and set a 10% tolerance for quadrupole 
effects. 4 



We begin with the sample of 28 early-type gravitational lens galaxies and bulges analyzed 
by Rusin et al. (2003). Note that this includes Q2237+050 (Huchra et al. 1985), which is lensed 
entirely by the bulge of a spiral galaxy. Twenty-two of the systems have measured lens redshifts 
(zd), and 18 of these are in combination with a measured source redshift (z s ). 5 Previously, for 
systems with no zj, the lens redshift had been estimated under the requirement that the galaxy fall 
on the fundamental plane (Kochanek et al. 2000). To make use of this technique, one must first 



4 We note in passing that when using steep radial profiles to fit fake lenses produced from shallower profiles, huge 
external shears or extremely flattened deflectors may be required. These solutions, however, can be rejected on the 
grounds of physical implausibility. First, the mass distributions of lens galaxies would have to be significantly more 
flattened than the light distributions, which typically have moderate axial ratios of ~ 0.7 — 0.8; (e.g., Keeton et al. 
1997). Second, such large external shears would have to be directly traced to nearby, cluster-scale perturbers, which 
are not observed in the relatively low-density environments of lens galaxies (Keeton, Christlein & Zabludoff 2000). 



5 Here we include a recent measurement of z s — 1.17 (Treu & Koopmans 2003) for MG1549+3047 (Lehar et al. 
1993). 
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estimate the stellar velocity dispersion, which necessitates an explicit assumption of a galaxy mass 
profile. Since this paper is concerned with constraining mass models, it is prudent to exclude the 
six lenses without spectroscopic lens redshifts. The number of lenses used in the following analysis 
is therefore 22. The systems are listed in Table 1. 

3.2. Geometric Properties 

The important geometric properties of each lens system are derived from the radii of the 
lensed images. Image positions with respect to the lens galaxy are determined from the HST data, 
based on the fitting methods outlined by Lehar et al. (2000). 6 The uncertainties on these positions 
are typically < 5 milliarcseconds, which is negligible compared to the radii of even the closest 
lensed images in doubles (~ 0.1 arcsec). For quads and rings, the angular Einstein radius TEin is 
determined by fitting an SIS plus external shear model to the lens data, and then converted to the 
physical radius Reiu = rEm-Drf- The Einstein radius can typically be derived to ~1% precision, so 
its uncertainty is also negligible. The geometric properties of the lens sample are listed in Table 1. 

Interpreting radii in terms of masses requires the critical surface density S cr , which depends on 
source and galaxy redshifts through the angular diameter distances. Recall that each lens galaxy 
in our sample has a measured redshift. If the source redshift has also been measured, then there 
is no uncertainty in S cr , at least within the context of our standard cosmological assumptions. If 
the source redshift is not known, we derive the uncertainty in S cr (actually log S cr ) with Monte 
Carlo techniques. Following Rusin et al. (2003), we draw 10000 values from a Gaussian distribution 
z s = 2.0± 1.0 for the source redshift, keeping only trials with 1 < z s < 5 because virtually all known 
lensed sources lie within this range. We adopt the z s = 2.0 value of log S cr as the median, and the 
rms scatter around this value as the uncertainty in logS cr . 

3.3. Photometric Properties 

The important photometric properties of each galaxy are the intermediate axis effective radius 
and luminosity. These quantities are derived from HST data obtained with the WFPC2 and 
NICMOS cameras. The fitting methods are detailed by Lehar et al. (2000), Kochanek et al. (2000) 
and Rusin et al. (2003). The optical effective radius (r e ; or R e = r e D<i in physical units) is 
determined by fitting a de Vaucouleurs profile to the galaxy in the filter with the highest signal-to- 



6 Because the two images of the radio lens B2319+051 (Rusin et al. 2001) have not been detected in the optical, a 
direct determination of their positions with respect to the galaxy is currently impossible. However, the structure in 
the lensed jets allows the galaxy position to be determined from mass modeling. Rusin et al. (2001) demonstrate this 
in the context of an isothermal mass profile, but our calculations show that the recovered galaxy position is nearly 
independent of the assumed model. We therefore include this system in our analysis, using the estimated Ri and R2 
from that paper. 
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noise ratio. This quantity is then held fixed to determine the mean surface brightness within the 
effective radius (n e ,Y ) in ah filters Y. The total magnitude is then my = fJ. e ,Y — 5 log r e — 2.5 log 2tt. 
The photometric properties of our sample are tabulated in Rusin et al. (2003). 

Complications arise from the need to convert the observed galaxy magnitudes to some standard 
luminosity scale that accounts for if-corrections and luminosity evolution. As we shall see in §4, 
only the magnitude offset M — M* = — 2.5 log(L/L*) is required for our analysis, since we can 
parameterize the homology model in terms of the present-day (z = 0) stellar mass-to-light ratio for 
an L* galaxy. Assuming that M — M* does not evolve (i.e., that the evolution rate is independent 
of luminosity; the same assumption was made in Rusin et al. 2003 and all other analyses using 
the FP to study galaxy evolution; e.g., van Dokkum et al. 1998, 2001; Treu et al. 2001, 2002), we 
can estimate this quantity from the observed magnitudes, without converting to rest frame bands. 
Specifically, we compare the observed magnitudes (m & Sj y) to models (m mo d,y) for the apparent 
magnitudes of an galaxy at the lens redshift, and then calculate the mean offset: 



(M-M,)= ^— y y^ru^Y) (3) 

l/(6m obSiY ) 2 

We compute the model magnitudes by convolving spectral energy distributions from the GISSEL96 
version of the Bruzual & Chariot (1993) spectral evolution models with transmission curves for HST 
filters (available from the technical archives of STScI, with zero-points from Holtzman et al. 1995). 
All models are normalized to a fixed present-day characteristic magnitude of M* = —19.9 + 51og/i 
for early- type galaxies in the B band (e.g., Madgwick et al. 2002). We take as a fiducial model 
an instantaneous starburst at Zf = 3 with solar metallicity Z = Z & and a Salpeter (1955) initial 
mass function (IMF). The assumed metallicity is consistent with observations of early- type field 
galaxies at 0.3 < z < 0.9 (Ferreras, Chariot & Silk 1999), and the assumed formation redshift 
is consistent with a number of analyses which strongly favor old stellar populations (mean star 
formation redshift (zf) > 2; e.g., Bernardi et al. 1998; Schade et al. 1999; Kochanek et al. 2000; 
van Dokkum et al. 2001; Im et al. 2002; van de Ven, van Dokkum & Franx 2002; Rusin et al. 2003). 
Considering a broad range of stellar models (1.5 < Zf < 5, 0.4 < Z/Zq < 2.5) yields an rms scatter 
in {M — M*) of 0.1-0.2 mag for most lenses. This uncertainty is much larger than the measurement 
errors, since all of the 22 galaxies have excellent photometry. We therefore adopt a luminosity error 
of 5{M - M*) = 0.20, or <51og(L/L*) = 0.08, for all lens galaxies. 

Estimating L/L* requires some understanding of luminosity evolution, which is implicitly 
included in the spectral models. For example, when normalized to a fixed L* at z = 0, a faster 
evolution rate (lower (zf)) would predict higher values of L*(z), and hence we would derive lower 
values of L/L* for each lens galaxy. If the true mean star formation redshift is (zf) < 1.5, then 
our fiducial model of (zf) = 3 could greatly overestimate L/L*, especially for high redshift lenses 
(zd > 0.7). However, most studies of early-type field galaxies are inconsistent with the rapid 
luminosity evolution implied by (zf) < 1.5. We note in passing that if the evolution rate is 
explicitly fit by our mass models (§4), the scatter is minimized by (zf) ~ 3, and increases for 
later formation redshifts. Unfortunately, the resulting constraints on the evolution rate are very 
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poor, as a number of other parameters are being fit. Since this paper focuses on the structure of 
lens galaxies, we believe it is best to simply assume a reasonable mean star formation redshift of 
(zf) = 3. We reiterate that our estimated uncertainties on L/L* are sufficient to account for the 
spread in evolution rates over the favored range of {zf) > 1.5. 

Finally, note that many of the quantities entering our analysis are correlated. For example, the 
total magnitude and the fraction of the light enclosed by the aperture radius R both correlate with 
the effective radius R e . Neither effect is very significant, because 5logR e is typically quite small 
(see Rusin et al. 2003), but for completeness we take these correlations into account via Monte 
Carlo calculations. We draw 10000 trials from a Gaussian distribution representing the effective 
radius (mean logi? e , width 5logR e ), compute the derived quantities, and store them for later use. 
We use the 10000 sets of correlated parameters, plus the independent (uncorrelated) luminosity 
error of <51og(L/L*) = 0.08, to calculate the scatter in the quantities entering the fitting functions 
given below. 

4. Self-similar Models for Galaxy Mass Distributions 

The aperture mass-radius relation (eq. 1 or 2) for a single lens tells us nothing about the mass 
profile. However, recall that the image splitting scale depends on the redshifts of both the lens 
galaxy (z^) and lensed source (z s ) through the angular diameter distances in S cr . For example, 
consider a spherical galaxy with a power-law mass density p(r) = po( r / r o)~ n i and a corresponding 
surface density S(i?) = T,q(R/ Rq) 1 ~ n . If a source sits directly behind this galaxy, it will be lensed 
into aring of Einstein radius Reiu = [2£o/X C r(3 — n)] 1 ^ n ~ 1 ^ Ro. We could therefore imagine placing 
a fixed mass distribution at various Zd, where it will lens sources at various z s . The resulting set 
of Einstein radii Reiu and aperture masses Meiu would then allow us to trace out the radial mass 
profile of that galaxy. 

Nature offers us an ensemble of gravitational lens systems with a wide range of lens (0 < < 1) 
and source (1 < z s < 5) redshifts, but the galaxy properties are far from uniform. The aperture 
mass-radius relations can, however, determine the typical (or "mean") mass profile of the lens galaxy 
population if all galaxies can be placed on a common mass scale. This is most easily accomplished by 
postulating some self-similar (homology) model to relate the mass and light distributions of early- 
type galaxies. A given model allows projected masses M{R) to be predicted for each lens using the 
photometric parameters, which can then be compared to the mass-radius relations determined from 
the lensing geometries. The large range of image radii in our sample, spanning 0.2 < R/R e < 7, 
gives us a long baseline for mapping the profile. We note, however, that the distribution is not 
uniformly populated, with most of the quads and rings probing the range 1 < REin/Re i$ 4. 
Doubles help fill in the radial coverage: the inner images typically span 0.4 < R\/R e < 1, while 
the outer images span 2 < R2/R e ^ 4. 

We consider a physically-motivated two component (luminous plus dark matter) model for the 
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mass distribution. Because of the limited sample of lenses on which to test the model, we keep 
the number of parameters to a minimum. The component which traces the luminous matter is 
modeled by two parameters: the mean present-day (z = 0) stellar mass-to-light ratio Y* (in the 
rest frame B band) for an L* galaxy, and an exponent (x) describing its dependence on the galaxy 
luminosity L: T oc L x . Under the assumption of homology, the fundamental plane implies that 
x ~ 0.3 (e.g., Jorgensen, Franx & Kjaergaard 1996). For each galaxy, the luminous surface mass 
density is represented by a de Vaucouleurs profile, and the associated mass enclosed by an aperture 
R is i 

M lum (R) = T*L* (J^j g (J^ , (4) 

where g(R/R e ) is the projected fraction of the luminosity inside R. Recall that for the de Vau- 
couleurs profile, g(0) = 0, g(l) = 0.5 and g(oo) = 1. Note that our use of the intermediate axis 
effective radius, which is the geometric mean of the major and minor axes, allows us to accurately 
determine the fraction of the luminosity enclosed by a circular aperture, even when the surface 
brightness distribution is elliptical. 

We model the dark matter as a power-law mass distribution. N-body simulations predict halos 
in which the mass density follows a shallow power-law slope at small radius, and turns over to a 
steeper slope beyond some break radius (e.g., Navarro, Frenk & White 1997, hereafter NFW; Moore 
et al. 1999b). However, because the lensed image radii probe only the inner several kiloparsecs of the 
mass distribution, we believe that a power-law approximation for the dark matter halo is reasonable 
on this scale. We describe this component by two parameters: the projected mass fraction (f cd m) 
within 2R e , which is roughly the median scale of the Einstein radii of our lens sample, and the 
logarithmic slope of the mass density profile (n, where p oc r~ n )J For each galaxy, the dark matter 
enclosed by an aperture R is 

M cdm (R) = M lum {2R e ) - (^-T n = T,L, ( f ) ^ g(2) - f^-V"" , (5) 

L-Jcdm \^tie) \ L *J 1 — Jcdm J 

where g(2) = 0.69. The total model-predicted mass inside R is then 

M mod (R) = M lum (R) + M cdm {R) . (6) 

Note that the parameters Y* and x set the normalization of the mass profile, while the dark matter 
parameters f c dm an d n modulate its shape. The homology scheme effectively places all galaxies on 
a common mass scale by normalizing the aperture masses M(R) by (L/L ! „) 1+:c , and the aperture 
radii R by R e . If a strict homology holds, then all galaxy aperture masses and radii should reside 
on a curve in the space of M{R)/{L/L lf ) l+X versus R/R e . 



7 The parameter n should be thought of as a "local" density slope. In reality, our fits use a surface density 
S oc R ~ n . This avoids the divergences encountered in converting between volume and surface densities for n < 1. 
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Modeling, dynamical and statistical studies of gravitational lens galaxies often make use of a 
single-component mass distribution with a scale-free radial profile. While simplistic, this model al- 
lows one to explore much of the lensing phenomenology related to mass concentration, and provides 
a standard for comparing profile constraints from different lens systems. It is therefore worthwhile 
to consider the scale-free limit of our homology model. By taking Y* — ► and f cdm — > 1 (such that 
T*/cdm/(l — fcdm) — ► const.), we obtain a pure power-law mass distribution: 



M pl (R) 



l+x 



R 

2R f 



3—n 



(7) 



where Mq is the projected mass inside 2R e for an L* galaxy. 



We will constrain the parameters of the homology model by the set of mass-radius relations 
from the lens sample. The optimization is performed using the goodness-of-fit function 



X 



2,2 
Xqr ' Xd " 



(8) 



In quads and rings, the lensing geometry constrains the mass enclosed by the Einstein radius (eq. 1) 
These systems are evaluated using 



Xqr 



E 



logM mod (R Einji ) - log[T, crti TrR 2 Ein i ] 

&scale&i 



(9) 



where 5 sca i e and Si are defined below. In doubles, the lensing geometry constrains a combination 
of two radii and projected masses (eq. 2). These systems are evaluated using 



x 2 d 



E 



\og[M mod {R hi )/R hi + M mod (R 2ii )/R 2:i ] - log[S cr ,i7r(fli,i + R 2 ,i)} 

&scale&i 



1 2 



(10) 



The logarithmic uncertainty <5j on each data point is derived using the Monte Carlo methods 
outlined in §3, but is well approximated as 5 2 ~ (1 + x) 2 (5 log L/L*) 2 + (5 log S cr ,i) 2 + ^ where 
5j is the additional 10% tolerance for quadrupole-related smearing of the mass-radius relation in 
doubles. We set 5 1 = for quads and rings, and 51ogS cr) j = for systems with a measured z s . 
Finally, we note that the above procedures have been tested using Monte Carlo simulations, and 
we find that the input mass profiles can be successfully recovered. 

Following optimization, we uniformly rescale the estimated errors by setting S sca i e so that 
the best-fit model has \ 2 = ^dof, the number of degrees of freedom. This does not alter the 
optimized parameters, but does allow us to relate the uncertainties to the observed scatter in 
the homology model. Moreover, it preserves the relative weighting among the data points, which 
naturally gives more weight to quads and rings, and those systems with a measured z s . Rescaling is 
particularly important because our model is undoubtedly a simplistic representation of galaxy mass 
distributions. Hence, the \ 2 is likely to be dominated by unmodeled complexity (i.e., deviations 
from self-similarity) in the galaxy population, rather than by observational errors. 
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5. Analysis and Results 

We first simultaneously constrain all four parameters (T#, x, fcdm-, n) in our homology model. 
The best-fit model has x 2 = 44.1 for N DO F = 18 (with 5 sca i e = 1). The rms scatter (0.15 in logM) 
is significantly larger than can be accounted for by the assumed observational uncertainties. Given 
our estimates of the measurement errors, the intrinsic scatter is roughly 30% in mass. Fig. 1 shows 
the pairwise parameter constraints. The contours correspond to A% 2 = 1, 2.30, 4 and 6.17. Unless 
otherwise noted, these and all subsequent Ax 2 values involve the rescaled x 2 > m which 5 sca i e is set 
so that the best-fit model has \ 2 = ^dof- 

There are well-defined regions allowed by the data, but substantial degeneracy between param- 
eters. This is expected because the mass-radius relations constrain the combined (luminous plus 
dark matter) radial mass distribution. For example, we require a projected mass of ~ 3 x 1O 11 M 
inside 2R e for the typical lens galaxy. This can be achieved by placing more mass in the lumi- 
nous component, which requires a larger Y*, or more mass in the dark matter component, which 
requires a smaller Y*. This degeneracy is reflected in the Y* — f c dm panel of Fig. 1. The mass-radius 
relations also prefer some overall shape, or concentration, for the combined profile. To achieve a 
given mass concentration, we can place more of the mass in the stellar component if the CDM slope 
is shallower, and less if the CDM slope is steeper. These degeneracies (Y* — n, f c ^ m — n) are also 
clearly seen in Fig. 1. 

The slope x, which describes the increase in stellar mass-to-light ratio with luminosity, is 
uncorrelated with the other three parameters (Fig. 1). Optimizing over the other parameters, 
we find x = 0.14^°;^ at 68% confidence (Ax 2 < 1), and -0.08 < x < 0.49 at 95% confidence 
(Ax 2 < 4). This is the first constraint on x from strong lensing. 

Our analysis robustly detects the presence of dark matter in the form of a mass component 
that is more spatially extended than the light. Dark matter contributes several tens of percent 
to the projected mass inside 2R e . This result can be quantified in a variety of ways, for different 
assumptions: 

• A model in which mass traces light (f c dm = 0) has Ax 2 = 10.7 with respect to the overall 
best fit, and is therefore rejected at > 99% confidence. It is too centrally concentrated to 
account for the ensemble of aperture mass-radius relations. We note in passing that such a 
model would require Y* ~ 11 Yq in the rest frame B band, only slightly higher than the value 
of Y* = (7.8 ± 2.7)Y determined by Gerhard et al. (2001) from dynamical analyses of the 
central regions of local (nearly L*) early- type galaxies. 

• Optimizing over the other three parameters, we find f c dm > 0.36 at 68% confidence, and 
fcdm > 0.08 at 95% confidence. 

• For a shallow CDM density slope of n = 1, as suggested by the NFW profile, the best-fit 
model has f c d m = 0.22, and Ax 2 = 2.4 with respect to the best overall model. The allowed 
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range is 0.12 < f cdm < 0.32 at 68% confidence, and 0.08 < f cdm < 0.45 at 95% confidence. 

• For a steeper CDM density slope of n = 1.5, as suggested by the Moore et al. (1999b) profile, 
the best-fit model has f c dm = 0.43, and A\ 2 = 0.8 with respect to the best overall model. 
The allowed range is 0.27 < f cdm < 0.57 at 68% confidence, and 0.14 < f cdm < 0.71 at 95% 
confidence. 

Because there are significant degeneracies between parameters, we cannot place a tight con- 
straint on the slope of the dark matter component: 1.44 < n < 2.20 at 68% confidence, with an 
upper bound of n < 2.33 at 95% confidence. However, the favored models have dark and luminous 
components which sum to produce very similar mass profiles over the radial range probed by the 
lensed images. This is illustrated in Fig. 2, which shows the projected mass M{R) versus R for all 
models with Ax 2 < 2.30 in the fcdm — n plane. We see that the mass profile closely tracks a pure 
isothermal (M oc R) model. 

In the scale-free limit (T* — > and fcdm — * 1) see §4), we find n = 2.07±0.13 at 68% confidence, 
consistent with isothermal or slightly steeper profiles. The 95% range is 1.83 < n < 2.33. The 
best-fit scale-free model has A% 2 = 0.1 relative to the best overall model, so it provides a good 
approximation to the mass distributions in lens galaxies. 

Fig. 3 demonstrates how the ensemble of mass-radius relations traces the radial mass profile. 
The solid line shows the scaled aperture mass M{R)/{L/L*) 1+X as a function of the scaled aperture 
radius R/R e , predicted by the best-fit power-law model. Quad and ring systems are plotted at a 
single radius Re™, with enclosed mass X] cr 7r R^in' Doubles are more difficult to depict, as they 
constrain a combination of two masses and radii. Note, however, that any mass profile can be 
normalized such that the mass-radius relation for a specific double is satisfied exactly. The resulting 
values of M{R\) and M{R^) i n the context of the best-fit profile are plotted for each double in 
Fig. 3. While the locations of the points obviously depend on the assumed model, the best-fit 
profile will minimize their scatter about the model curve. In this way we see that doubles also help 
map the mass profile. 

The radial coverage offered by the current lens sample is not uniform (Fig. 3). Quad and ring 
systems, which provide the strongest mass constraints, mostly lie in the range 1 < REin/Re ^5 4. 
There are two notable outliers: Q2237+030 (Huchra et al. 1985), which probes the smallest radial 
scale {REin/Re — 0.2), and MG2016+112 (Lawrence et al. 1984), which probes the largest radial 
scale {REin/Re — 7). Because these lenses constrain the extremes of the density profile, it is 
important to address how much they affect our fits. We find that dropping either of the systems does 
not alter the primary conclusions that the total mass profile is nearly isothermal, or that x ~ 0.14 is 
favored. Some of the other constraints are slightly affected, and this can be understood based on the 
unique characteristics of each lens. First, MG2016+112 probes a scale much larger than the effective 
radius, and its aperture mass will be dominated by dark matter in most models. This system is 
therefore important for constraining the properties of the CDM component. Removing it weakens 
the exclusion of models in which mass traces light to about 95% confidence, and significantly 
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weakens the upper limit on n in the limit fcdm ~^ 0. Profile constraints in the scale-free limit 
(fcdm ~^ 1)) however, are not significantly changed. Second, Q2237+030 probes a scale much 
smaller than the effective radius, and its aperture mass will be dominated by luminous matter 
in most models. Consequently, this system contributes minimally to constraints on the CDM 
component. Dropping it has little effect on any of the model parameters, even those related to the 
stellar mass. 

Finally, it is interesting to dissect the distribution of scatter about the best-fit homology 
model. Despite comprising only 7 of the 22 systems, doubles account for 56% of the x 2 - The 
additional scatter for doubles is greater than we would expect from the quadrupole smearing, and 
may indicate some unmodeled effect which we do not presently understand. Dropping all doubles 
would significantly improve the fit (x 2 = 17.1 for N DOF = 11, 

prior to rescaling) and tighten the 
resulting parameter constraints, but would not alter the major conclusions. For example, a model 
in which mass traces light would be rejected more strongly (Ax 2 = 19.6), and constraints on the 
density slope in the scale-free limit would be n = 2.04 ± 0.09 (68% C.L.). We also note that a 
pair of lenses are particularly strong contributors to the scatter. The two-image system Q0142-100 
(Surdej et al. 1987) is the farthest outlier, just as it is in our FP analysis (Rusin et al. 2003), and 
accounts for almost 1/4 of the x 2 - Because the galaxy appears to be much brighter than expected, 
it greatly over-predicts aperture masses. The ring system MG1131+0456 (Hewitt et al. 1988) is 
the second largest contributor to the x 2 , accounting for 1/7 of its value, despite being assigned a 
rather large fit tolerance due to its estimated source redshift. The lens galaxy can be brought into 
better agreement with the other galaxies if the lensed source is at a substantially higher redshift 
than the assumed median of z s = 2. While our sample is too small to reasonably remove these 
two outliers, doing so would improve the fit (x 2 = 28.3 for N DOF = 16, prior to rescaling) but not 
affect the major conclusions. For example, a model in which mass traces light would be rejected 
somewhat more strongly (Ax 2 = 11.1), and constraints on the density slope in the scale-free limit 
would be n = 2.08 ± 0.10 (68% C.L.). 



6. Applications 

6.1. Time Delays and the Hubble Constant 

The time delays between gravitationally lensed images depend on a combination of the Hubble 
constant and the mass distribution of the lensing galaxy. The Hubble constant inferred from a 
measured time delay can be approximated as 

H = A(1-(k))+B{k)(ti-1) + C , (11) 

where (k) = (£)/£ cr is the mean scaled surface mass density in the annulus defined by the radii 
of the lensed images, r] is the logarithmic slope of the density (k oc R 1 ~ ri ) within the annulus, 
and the constants A and B incorporate the lensing geometry, the redshifts, and the measured 
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delay (Kochanek 2002). The A term is the most important, and the B term just contributes a 
small (~10%) correction. Kochanek (2002, 2003) investigates "simple" time delay lenses which are 
dominated by a single galaxy with a precisely determined centroid relative to the images. Two 
limiting cases of the mass distribution are considered: models in which mass traces light yield 
h ~ 0.7, while isothermal models yield h ~ 0.5. A number of recent results from a broader 
(perhaps less well understood) sample of lenses fall within the range 0.45 < h < 0.65 (Burud et al. 
2002a, 2002b; Fassnacht et al. 2002; Treu & Koopmans 2002b; Winn et al. 2002) if isothermality 
is assumed. In fact, the only lens to yield h ~ 0.7 for an isothermal model is B0218+357 (Biggs 
et al. 1999), but the result is currently dominated by systematics related to the poorly constrained 
galaxy position. One can therefore make the case that estimates of the Hubble constant from 
strong lensing are in broad agreement, but are systematically lower than the values favored by 
both the HST Key Project (Freedman et al. 2001) and the recent WMAP results (Spergel et 
al. 2003). Because Kochanek (2002, 2003) shows that nearly constant mass-to-light ratio models 
are necessary to reconcile time delay measurements with h ~ 0.7, this would seem to imply that 
lens galaxies have no extended dark matter component, in contradiction to virtually every other 
observational test. 

We can now explore this intriguing problem in the context of our constraints on self-similar 
galaxy models. Of the five lenses considered by Kochanek (2002), three are drawn from our sample. 
The two exceptions are B1600+434 (Jackson et al. 1995), which has a late-type lens galaxy (Koop- 
mans, de Bruyn & Jackson 1998) and is therefore not described by our models, and RXJ091 1+0551 
(Bade et al. 1997), which has a significant cluster perturbation complicating interpretation of the 
image radii. Because PG1115+080 (Weymann et al. 1980) has been explored in detail by Treu 
&: Koopmans (2002b), we limit our demonstration to a pair of two-image lenses, SBS1520+530 
(Chavushyan et al. 1997) and HE2149-2745 (Wisotzki et al. 1996). 

Because the favored homology models produce very similar mass profiles, (k), ij and the implied 
Hq do not vary much along the f c dm — n degeneracy stripe. These three quantities are plotted 
for SBS1520+530 and HE2149-2745 in Figs. 4 and 5, respectively, by normalizing the models to 
reproduce the mass-radius relations exactly. The values of Hq use the coefficients A, B and C from 
Kochanek (2002). The uncertainties in these coefficients are propagated to Hq using Monte Carlo 
procedures. Assuming that each model {f c dm,ii n j) has a likelihood pij oc exp(— A%? /2), constraints 
on Hq can be determined from a Bayesian analysis. These constraints are shown in panel (d) of 
Figs. 4 and 5. For SBS1520+530, h = 0.58 ± 0.08 (68% C.L.). For HE2149-2745, h = 0.55 ± 0.10 
(68% C.L.). The values of h come out slightly higher than those obtained by Kochanek (2002) for 
an isothermal profile, because our favored profiles are slightly steeper than n = 2. The Hubble 
constants are still lower than those favored by the HST Key Project (h = 0.72 ± 0.08; Freedman et 
al. 2001) and WMAP data (h = 0.72 ±0.05; Spergel et al. 2003), but the discrepancy is only about 
la in our analysis. 



6.2. Flux Ratio Anomalies and Microlensing 



Simple, smooth models for the galaxy mass distribution have great difficulty reproducing the 
observed flux ratios in four-image lenses (e.g., Dalai & Kochanek 2002; Metcalf & Zhao 2002; Chiba 
2002; Keeton et al. 2003). Moreover, flux ratio measurements in both the radio and optical appear 
to violate fundamental symmetry arguments which predict nearly equal magnifications for merging 
image pairs. Such flux ratio "anomalies" are almost certainly due to small-scale structure in the 
gravitational potential (Kochanek & Dalai 2003), which can significantly perturb magnifications 
(e.g., Mao & Schneider 1998) while leaving the image positions virtually unaltered. There is a debate 
regarding the identification of the perturbers. CDM simulations (Moore et al. 1999a) predict sig- 
nificant substructure (with mass scale M > 1O 6 M ) in galaxy halos, and Dalai & Kochanek (2002) 
demonstrate that the substructure mass fraction required to account for radio flux ratio anoma- 
lies is in line with iV-body results. Stars are another source of substructure in the gravitational 
potential (Witt et al. 1995), and microlensing- induced time variability has been unambiguously 
detected in the light curves of Q2237+030 (Wozniak et al. 2000). Schechter & Wambsganss (2002) 
have pointed out that there is a tendency for saddle point images to be demagnified compared to 
our expectations for smooth lens models (see also Kochanek & Dalai 2003). They demonstrate 
that stellar microlensing produces this effect only if the stars account for a relatively small fraction 
(0.15 < Kium/K ^$ 0.30) of the surface mass density at the image radius. 8 If Ki U m/ K is too high, 
then the asymmetric effects on saddles and local minima disappear; if Ki um /K is too low, then the 
anomalies are too rare. We can now test whether our two-component mass models fall within the 
favored range. 

In the context of our homology model, the value of Ki um / k at a given radius R/R e depends only 
on the dark matter abundance parameter f c dm and density slope n. Fig. 6a plots ni um /n at 2R e , 
a typical radial distance for lensed images. Since we require models with an extended dark matter 
component, lower values of Ki um /n are favored at 2R e due to the fact that the surface brightness of 
the galaxy has already decreased substantially at this radius, while the dark matter density has not. 
Assuming that each model (f c dm,i, n j) has a likelihood pij oc exp(— Ax?-/2), constraints on Ki urn /n 
can be determined from a Bayesian analysis. These constraints are shown in Fig. 6b. At R e , just 
about any value of Ki um /n is permitted. At 2R e , however, the results are quite restrictive, favoring 
% ra //t < 0.31 (68% C.L.). Hence, our models are consistent with the values of Ki um /K needed to 
produce the relative demagnification of saddle point images observed in flux ratio anomalies. 



8 To ward off possible confusion of terminology, we note that while f cdm is a global parameter describing the 
projected dark matter mass fraction inside 2R e , Kx um j 'k is a local parameter describing the fraction of surface mass 
density in the form of stars at the radius of a lensed image. 



7. Discussion 



We have constrained the typical mass profile of early-type lens galaxies using an ensemble of 
aperture mass-radius relations from 22 gravitational lenses. The different galaxies are combined 
using a self-similar mass model consisting of four parameters: the present-day (z = 0) stellar 
mass-to-light ratio in the B band for an L* galaxy (T*), its dependence on galaxy luminosity (x, 
where T oc L x ), the projected mass fraction of dark matter inside 2 effective radii (f c dm), an d 
the logarithmic density slope of the CDM (n, where p oc r~ n ). Despite significant parameter 
degeneracies, the favored models have dark and luminous components which sum to produce very 
similar total mass profiles over the radial range probed by the lensed images. In the scale-free limit 
we find n = 2.07 ± 0.13 (68% C.L.), consistent with isothermal or slightly steeper profiles. Two- 
component models imply that the dark matter mass fraction is not very high: 0.12 < f c d m < 0.32 
for an n = 1 halo (NFW 1997), or 0.27 < f cdm < 0.57 for a steeper n = 1.5 halo (Moore et al. 
1999b), both at 68% confidence. Even so, the need for an extended dark matter component in 
the inner regions of the galaxies is significant, as a model in which mass traces light is ruled out 
at > 99% confidence. These are the best statistical constraints yet obtained on the typical radial 
mass profile in gravitational lens galaxies, and they add to evidence that the mass distribution is 
nearly isothermal on the scale of a few effective radii. 

One might expect lenses to be a biased sample of early-type galaxies, because deflectors with 
more centrally-concentrated profiles have larger lensing cross sections per unit mass. Consequently, 
if the galaxy population exhibits a range of profiles, then lensing would tend to select galaxies 
with steeper mass distributions. Note, however, that profile constraints for the lens galaxy pop- 
ulation are consistent with the nearly isothermal profiles favored by dynamical and X-ray studies 
of local ellipticals, which are not selected on the basis of mass concentration. This would appear 
to strengthen the argument that early-type galaxies do not exhibit much structural variety, and 
therefore, that gravitational lenses are typical members of this galaxy population. Of course, even if 
lenses were a biased galaxy sample, it would have no effect on applications such as Hubble constant 
determination, which only requires an understanding of the mass profiles in lens galaxies. 

Our results have implications for the relationship between dark and luminous matter in early- 
type galaxies, and can be compared with recent studies involving the stellar dynamics of local 
galaxies (Gerhard et al. 2001) and lenses at intermediate redshift (KT), the number and size 
distribution of lenses (Keeton 2001), and the existence of a tight fundamental plane (Bordello, 
Salucci k, Danese 2003). While the use of different mass models and priors makes a detailed 
comparison difficult, we can nonetheless survey the basic results for consistency. 

First, consider evidence for the existence and abundance of dark matter within a few optical 
radii. Our results rule out a galaxy population in which mass traces light at > 99% confidence, 
and imply a small but significant dark matter mass fraction inside 2R e . KT likewise reject con- 
stant mass-to- light ratio models based on the measured velocity dispersions of lens galaxies. In 
comparison, Gerhard et al. (2001) claim 10-40% dark matter within the volume defined by R e , but 
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caution that some elliptical galaxies show no evidence for dark matter at this scale. Coming from 
the other side, Keeton (2001) uses a variety of lensing statistics to place an upper limit of 40% dark 
matter within the volume defined by 2R e . Bordello et al. (2003) also suggest that a cuspy dark 
matter halo can contribute little to the mass inside a few effective radii, if early-type galaxies are to 
occupy a tight fundamental plane. We can compare these measurements to the predictions of the 
CDM model. Assuming an NFW halo profile, a cooled baryon fraction of f coo i = £lb,cooi/^M — 0.02 
(from the local census of baryons by Fukugita, Hogan k Peebles 1998), and a typical initial halo 
concentration (c ~ 9; Bullock et al. 2001), the adiabatic compression models of Keeton (2001) 
predict f c d m — 0.6-0.8. The predicted value of fcdm can be reduced slightly by increasing the 
cooled baryon fraction, but even maximally efficient cooling (f CO oi = Qb/^M = 0.17; Spergel et al. 
2003) yields fd m — 0.4-0.5. Because the above suite of measurements seem to imply lower dark 
matter mass fractions in the central regions of early-type galaxies, they may support the evidence 
from late-type galaxies (e.g., McGaugh k de Blok 1998; de Blok et al. 2001; Salucci 2001; de Blok 
k Bosma 2002) that halos are less concentrated than predicted by numerical simulations of CDM 
(e.g., Moore et al. 1999b; Bullock et al. 2001). 

q i 4 +0.16 

Next consider the dependence of mass-to-light ratio on luminosity. We find that T oc L ' -° 12 
(68% C.L.), the first such constraint from strong lensing. While this result strictly applies to the 
stellar mass-to-light ratio, the self-similarity of our model means that it also applies to the total 
mass-to-light ratio inside any fixed fraction of the effective radius. The luminosity dependence 
is slightly shallower than the slope of 0.5 < x < 0.7 measured by Gerhard et al. (2001) in the 
B band using stellar dynamics. Each of these results is broadly consistent with the value of 
x ~ 0.30 ± 0.05 obtained by interpreting the -B-band FP under the assumption of structural 
homology (e.g., Jorgensen et al. 1996). Similar investigations have recently been performed at 
longer wavelengths. Bernardi et al. (2003) measure x = 0.14 ± 0.02 in the r* band using virial 
mass-to-light ratios, and derive a similar value from their FP slopes. In addition, Borriello et al. 
(2003) show that scatter in the local r-band FP is minimized by x ~ 0.2. Finally, it is interesting 
to note that the above estimates of x, determined at the scale of the optical radius, are similar 
to the relations between virial mass and luminosity determined from weak lensing and satellite 
dynamics: Guzik k Seljak (2002) and Prada et al. (2003) find M vir oc L L2 " L6 , while McKay et al. 
(2001, 2002) favor M v i r oc L. The broad consistency of these scaling laws suggests that the mass 
distribution in early-type galaxies is closely related to the light over many orders of magnitude in 
radius. 

We have applied our constraints on the mass distribution of early-type galaxies to a pair of 
interesting problems. First, we find that the measured time delays in SBS1520+530 and HE2149- 
2745 favor a Hubble constant of h ~ 0.55 — 0.60. This value is slightly higher than that derived by 
Kochanek (2002) using isothermal models, since we prefer mass profiles that are slightly steeper 
than n = 2. The derived Hq is still systematically lower than recent results (Freedman et al. 2001; 
Spergel et al. 2003), although the discrepancy is only about la in our analysis, as the statistical 
uncertainties in the mass profile are still significant. Second, we have estimated the fraction of 
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the surface mass density in the form of stars (ni um /K) at the radii of lensed images. For images 
at 2i? e , our models favor ni um /n < 0.31 (68% C.L.). This is consistent with the values at which 
stellar microlensing can reproduce the observed flux ratio anomalies, in particular the tendency for 
saddle-point images to be more strongly perturbed than minima (Schechter & Wambsganss 2002). 

Much more work is needed to fully understand the structure of early-type galaxies, and the 
implications for galaxy formation theories. Because of its unique advantages, gravitational lensing 
is certain to contribute significantly to this goal. Improved statistical tests to constrain the mass 
profiles of early-type galaxies require larger samples of lenses. Particularly important is the addition 
of new quad and ring systems which probe the mass distribution on small (REin < Re) and 
large (REin > 4i? e ) scales, and therefore fill in the tails of our radial coverage. As always, the 
key to turning recently discovered systems into useful astrophysical tools is the acquisition of 
high-quality HST photometry and ground-based spectroscopy. Expanded samples could allow us 
to replace our simplistic power-law model for the dark matter halo with more realistic profiles, 
and consider structural dependencies on luminosity or color. In this way we may be able to 
properly investigate sources of scatter in our homology model, and perhaps quantify or constrain 
the structural diversity of the early-type galaxy population. With regard to lensing determinations 
of the Hubble constant, it is vital to increase the number of lenses that have both direct profile 
constraints and well-determined time delays. Currently there is only one such system (PG1115+080; 
Schechter et al. 1997; Treu & Koopmans 2002b). 9 Programs are now underway to both measure 
more time delays, and obtain profile constraints on current systems by measuring stellar velocity 
dispersions or mapping extended emission from the lensed host galaxies. Such observations should 
greatly illuminate current discrepancies related to time delays and the Hubble constant, and, more 
generally, improve our understanding of the mass distribution in early-type galaxies. 

We thank Josh Winn and the anonymous referee for offering comments and suggestions which 
greatly improved this manuscript. We acknowledge the support of HST grants GO-7495, 7887, 8175, 
8804, and 9133. We acknowledge the support of the Smithsonian Institution. CSK is supported by 
NASA ATP Grant NAG5-9265. 



9 A number of profile constraints have been claimed for the time delay lens Q0957+561 (e.g., Grogin & Narayan 
1996), but Keeton et al. (2000) demonstrates that existing models fail to properly reproduce lensed extended emission 
from the quasar host galaxy. 
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Table 1. Structural Parameters for Lens Systems 



Lens 




z s 


Morph. 




Radii (") 


0047-2808 


0.49 


3.60 


Quad 




REin = 1-35 


Q0142-100 


0.49 


2.72 


Double 


Ri 


= 0.38, R 2 = 1.86 


MG0414+0534 


0.96 


2.64 


Quad 




REin = 1-19 


B0712+472 


0.41 


1.34 


Quad 




R Ein = 0-71 


HS0818+1227 


0.39 


3.12 


Double 


Ri 


= 0.61, R 2 = 2.22 


B1030+074 


0.60 


1.54 


Double 


Ri 


= 0.18, R 2 = 1.39 


HE1104-1805 


0.73 


2.32 


Double 


Ri 


= 1.10, R 2 = 2.09 


PG1115+080 


0.31 


1.72 


Quad 




REin = 1-15 


MG1131+0456 


0.84 


- 


Ring 




REin = 1-05 


HST14113+5211 


0.46 


2.81 


Quad 




R E in = 0.86 


HST14176+5226 


0.81 


3.40 


Quad 




R E in = 1-42 


B1422+231 


0.34 


3.62 


Quad 




REin = 0.78 


SBS1520+530 


0.72 


1.86 


Double 


Ri 


= 0.39, R 2 = 1.21 


MG1549+3047 


0.11 


1.17 


Ring 




REin = 115 


B1608+656 


0.63 


1.39 


Quad 




R E in = 1-14 


MG1654+1346 


0.25 


1.74 


Ring 




REin = 1-05 


B1938+666 


0.88 




Ring 




R E in = 0.50 


MG2016+112 


1.00 


3.27 


Quad 




REin = 1-63 


B2045+265 


0.87 




Quad 




R Ein = 1-14 


HE2149-2745 


0.50 


2.03 


Double 


Ri 


= 0.34, R 2 = 1.35 


Q2237+030 


0.04 


1.69 


Quad 




REin = 0.88 


B2319+051 


0.62 




Double 


Ri 


= 0.61, R 2 = 0.82 



Note. — Listed for each lens are the galaxy (zd) and source (z s ) 
redshifts, morphology, and relevant radii (Einstein radius R E i n for 
quads and rings; image radii R\ and R 2 for doubles). Uncertainties 
on these radii are negligible. The image radii for B2319+051 are 
the estimated values from Rusin et al. (2001). 
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CDM mass fraction (<2R e ) f cdm CDM density slope n 

Fig. 1. — Constraints on the self-similar model. We show pairwise constraints on the present-day 
(z = 0) stellar mass-to-light ratio in the B band for an L* galaxy (T*, in solar units), its logarithmic 
dependence on luminosity (x), the projected CDM mass fraction inside 2 effective radii {fcdm)-, an d 
the logarithmic CDM density slope (n). Solid contours represent Ax 2 = 2.30 and 6.17, the 68% 
and 95% confidence levels for two parameters. Dotted lines represent A\ 2 = 1 and 4, the 68% and 
95% confidence levels for one parameter. The errors have been rescaled so that the best-fit model 
has x 2 = N DO f- 



-26- 



CM 



cm 



CO 
CO 
D 

E 

0) 

O 
CD 



7? O 



O 
CL 



1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 

y 




s 


y 

/ 




/ 

/ 

y 




/ • ■ .jidKr 




V 

I 1 I 






0.1 



10 



R/R 



Fig. 2. — Mass profiles for allowed two-component homology models. Each model falls within the 
68% (Ax 2 < 2.30) confidence region of the f c d m — n plane (Fig. 1). Solid lines are the projected 
masses inside R/R e , where R e is the effective radius. Profiles are normalized to a fixed projected 
mass at R = 2R e . For comparison we show the de Vaucouleurs profile (dotted line), and an offset 
isothermal profile (dashed line). While the allowed models exhibit a wide range of dark matter 
abundances, they all have total mass profiles which are approximately isothermal over the radial 
range spanned by the lensed images. 
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Fig. 3. — The aperture mass profile in the power-law limit. We plot the scaled projected mass 
M(R)/(L/ ' L^) l+X as a function of the scaled radius R/R e . The solid line is the prediction of 
the best-fit power-law model with n = 2.07 and x = 0.14. Data points represent individual lens 
systems. Solid squares are lenses with measured source redshifts; open squares are lenses with 
estimated (z s = 2.0 ± 1.0) source redshifts. For quads and rings, a single point is plotted showing 
REin and mass Yj cr iiR 2 Ein . For doubles (big circles), points are plotted at each image radius [R\ 
and R2) and connected with a dotted line, and the masses are the M{R\) and M(i?2) which would 
satisfy the mass-radius relation exactly in the context of this profile (see §5). 
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CDM mass fraction (<2R g ) f cdm H Q (km/s/Mpc) 

Fig. 4. — Interpreting the time delays of SBS1520+530. The first three panels show various 
quantities as a function of CDM mass fraction fdm an d CDM density slope n. Solid contours 
indicate (a) the mean scaled surface density (k) inside the image annulus, (b) the effective density 
slope (rj) inside the image annulus, and (c) the derived Hubble constant Hq. The models are 
normalized to satisfy the mass-radius relation exactly. Note that the quantities exhibit the expected 
behavior: rj — > n for fcdm 1 and (n) = 0.5 for a scale-free isothermal model. Also, we see that an 
isothermal profile yields Hq ~ 50kms -1 Mpc -1 , as demonstrated by Kochanek (2002). The dotted 
contours show the 68% (A% 2 < 2.30) and 95% (A% 2 < 6.17) confidence regions from the homology 
model (Fig. 1). Plotted in panel (d) is the relative probability of Hq based on the information 
displayed in (c). A Hubble constant of Hq = (58 ± 8) kms" 1 Mpc" 1 is favored (68% C.L.). 
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Fig. 5. — Interpreting the time delays of HE2149-2745. The plots are analogous to those in Fig. 4. 
A Hubble constant of H = (55 ± 10) kins" 1 Mpc -1 is favored (68% C.L.). 
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Fig. 6. — The fraction of the surface mass density in the form of stars, (a) Plotted are contours 
of Kium/K at R = 2R e for models in the fdm — n plane. The dotted contours show the 68% 
(Ax 2 < 2.30) and 95% (Ax 2 < 6.17) confidence regions from the homology model (Fig. 1). (b) 
The relative probability of ni um /n based on the information displayed in (a). The solid line is the 
fraction at R = 2R e . The dashed line is the fraction at R = R e . 



